Electromagnetic radiation in a time-varying background medium 

Neil V. Budko 

Laboratory of Electromagnetic Research, Faculty of Electrical Engineering, 
Mathematics and Computer Science, Delft University of Technology, 
Mekelweg 4, 2628 CD Delft, The Netherlands* 
(Dated: June 2, 2009) 

Abstract 

Analytical solutions are presented for the electromagnetic radiation by an arbitrary pulsed source 
into a homogeneous time- varying background medium. In the constant-impedance case an explicit 
radiation formula is obtained for the synchronous permittivity and permeability described by any 
positive function of time. As might be expected, such a medium introduces significant spectral 
shifts and spatio-temporal modulation, which are analized here for the linear and exponential time- 
variations of the medium parameters. In the varying-impedance case the solution is obtained for 
the fourth-order polynomial time-dependence of the permittivity. In addition to the spectral shifts 
and modulation this spatially homogeneous medium scatters the field introducing causal echoes at 
the receiver location. 
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I. INTRODUCTION 



One of the most difficult and least understood topics in science is the interaction between 
dynamic processes with separate physical origins and causes. Yet, such processes are very 
common in nature and engineering. For example, the electromagnetic field interacts with 
a non-autonomous (time-varying and independently controlled) medium whenever we try 
to visualize or measure an ultra-fast thermal, chemical or mechanical transition. A non- 
automous medium is also a good approximation for the nonlinear medium with a strong 
pumping field and a weak probing signal. Potential applications include: novel imaging 
algorithms for the visualization of ultra-fast dynamic processes in non-autonomous physi- 
cal systems undergoing mechanical displacement, deformation, phase transition, chemical 
reaction, etc; detection of weak transient variations in the gravitational field; creation of 
ultra-short impulses and spatially localized light bullets; fast adaptive optical elements; 
broad-band extension of photonic and metamaterial crystals and associated effects, such as 
photonic band-gap waveguides and switches, perfect lenses, and invisibility cloaks; creation 
of congestion-free and transient-proof power, information, and transport networks. 

At present we are generally aware that the electromagnetic field evolution in a non- 
autonomous system cannot be reduced to a sequence of snapshots corresponding to the 
instantaneous values of the system parameters at various times. A plethora of very inter- 
esting and poorly understood transient effects is expected to be present. Some of those are: 
cross-modulation, velocity modulation, and parametric amplification of waves, e.g. dramatic 
pulse reshaping in attosecond optics, and exponential growth of disturbances in networks. 
So far, however, we do not have any general theoretical picture of these phenomena and 
do not understand connections between them, and therefore do not yet realize their full 
technological potential. Most of the practical modeling is numerical and does not allow for 
in-depth analysis. A few available, mainly one-dimensional, analytical models [1, 2, 3] are 
often quasi- stationary, perturbative, restrict the type of the time-variation of the system 
parameters and the source, and have little to say about the near-field effects. Practically 
nothing is known about the influence of non-autonomous systems on the coherence of the 
field. Three-dimensional studies [4, 5] were mostly concerned with the closed systems, such 
as electromagnetic cavities and waveguides with perfectly conducting walls. In general, 
the authors of non-perturbative studies agree that the time-variation of the medium at a 
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scale comparable to the time-variation of the source results in a significant spatio-temporal 
modulation of the field. Recently there has been some progress in the theory of media period- 
ically varying in time [6], where, as it turns out, one can apply the temporal analogue of the 
Bloch theorem. The analysis shows that such time- varying media introduce bandgaps in the 
wavenumber (momentum) domain, as opposed to the frequency-domain (energy) bandgaps 
produced by the spatially periodic media. 

To understand the spatio-temporal modulation of the time-domain source by the gen- 
eral time-varying medium one has to carry out a systematic theoretical analysis of the 
simplest non-autonomous electromagnetic system - a homogeneous isotropic time-varying 
background medium. This is the subject of the present paper, where we shall consider the 
realistic vectorial three-dimensional field and arbitrary time-varying sources. This study is 
intended to fill some of the gaps in our knowledge, e.g. about the near-field behaviour, po- 
larization effects, and a few-cycle pulsed sources. In addition, it provides valuable analytical 
soluitions, useful for the verification of numerical codes. The mathematical techniques ap- 
plied here are standard, albeit involved. First, using a straightforward change of variable we 
obtain a completely explicit analytical solution for the constant-impedance case, where the 
time-variation of the permittivity and permeability is synchronous, but otherwise arbitrary 
(positive function of time). The more difficult time- varying impedance case is treated using 
an adaptation of the method originally proposed by Shvartsburg which allows to derive an 
explicit solution for a polynomial time-dependence of the permittivity. 

II. THE CONSTANT IMPEDANCE CASE 

We begin with a relatively simple problem where the permittivity and pemeablity of the 
medium vary synchronously, i.e. 



so that the medium impedance remains constant at all times, e{t)/ix{t) = £0/^0 = const. 
Consider Maxwell's equations 



e(t) = e a(t), fx(t) = fx a(t), 



(1) 



V x H(x,t) + d t D(x,t) = - 

V x E(x,i) + d t B(x,t) = 0, 



J(x,t), 



(2) 
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with the corresponding constitutive relations 

D(x, t) = e a(*)E(x, t), B(x, t) = /i a(t)H(x, t). (3) 

Our goal is to obtain an explicit solution, called radiation formula, which gives the fields in 
terms of the given source current density J. Although, the procedure is pretty standard, it 
will be presented in full detail, if only to clarify the differences with the case of time- varying 
impedance discussed in the next section. First, we rewrite the Maxwell equations as 

-V x B(x,t) + / u a(t)<9 t D(x,t) = -fi a(t) J(x, t), 

(4) 

V x D(x,t) +e a(t)d t B(x,t) = 0. 
Now, we introduce a new variable with the physical dimension of time: 

to 

The original time t will now be considered an implicit function t(r) defined by the expression 
above. Taking into account that 

d dr d 1 d . . 

we rewrite the Maxwell equations in terms of variable r as 

-V x B(x, r) + // d r D(x, r) = -/x a(r) J(x, r), 

(7) 

V x D(x, r) + £ <9 T B(x, r) = 0, 

where the known functions a(r) and J(x, r) must be understood as a(t(r)) and J(x, t(r)). 
The function t(r) is the inverse mapping with respect to the substiution (5). This inverse 
mapping exists and is one-to-one, if a(t') > 0, t < t' < t. Maxwell's equations (7) are 
identical to the Maxwell equations in vacuum and have a well-known analytical solution. 
One usually starts by eliminating one of the unknowns, say B, and arriving at the following 
second-order vector wave equation for D: 

V x V x D(x, r) + ^D(x, r) = -\d T [a(r) J(x, r)] , (8) 

where Co = l/^/e fi . From (7) we derive the compatibility relation: 

r 

V • D(x, r) = - (V- ) J a(r') J(x, r') dr'. (9) 
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Rewriting (8) as 



V (V • D(x, r)) - AD(x, r) + -<9 T 2 D(x, r) = --d T [a(r) J(x, r)] , 
and applying (9), we arrive at 

r 

AD(x,r)-^D(x,r) = 4«9 T Ji(x,r)-V(V-) / J 1 (x,r')dr', 



(10) 



(11) 



where we have introduced 



Ji(x, r) = a(r)J(x, r), 



(12) 



Performing the Laplace transform with respect to r and the three-dimensional Fourier trans- 
form with respect to x we obtain the following algebraic equation in the (k, s)-domain 



-|k| 2 D(k, s) - 7oD(k, s) = -J^k, s) + k(k- )-Ji(k, s), 



c5 s 



where 7 = s^fe^jM) = s/co- This equation can be solved as 

1 



D(k, s) 



|k| 2 + 7o 2 



4ji(k, s) + k(k-)-Ji(k, s) 



(13) 



(14) 



Carrying out the inverse three-dimensional Fourier transform we obtain the (x, s)-domain 
solution as 



D(x,s) = - J ( ? (x-x / , S )^J 1 (x', S )dx / + V(V-) j g{ 



x-x , ,s)-Ji(x , ,s)dx', (15) 

s 



where the scalar Green's function is 



#(x, s) 



e -7o|x| 

47r|x| ' 



(16) 



Computing the spatial derivatives and grouping the terms according to their spatial decay 
factors we get 



D(x,s) 



x'ei 



+ 



x'ei 



+ 



'eR 3 



e -7o|x- 


-x'| 


47T X — 


x'| 3 


e -7o|x- 


-x'l 


Att X — 


x'| 2 


e -7o|x- 


x'l 


4ir x — 


x'l 



I] -Ji(x',s) dx' 

s 



'if 

Co 
S i 



(17) 
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where Q = 0(0- ), = (x — x')/|x — x'| and I is the Kronecker identity tensor. The inverse 
Laplace transform of expressions (15) and (17) gives the (x, r)-domain solution 



gk 3 



E(x,t) = 

+ / -[3Q-^|— — — -dx' ( ;i8) 

1 " lv Y "- c € a(T) 



'GM 3 



1 


47T X - 
1 


x'| 3 


47T X — 
1 


x'| 2 


4-7T X — 


x1 [ 


e(r)E( 


x,r) 




= r 



/" 1 j| 5 r [a(r)J(x',r)] T=rr , 

,/ 47r|x — x'| CQ£ o(' r ) 



x — X 



(19) 

Co 

To obtain the (x, t)-domain solution we now need to introduce the inverse r-function t(r), 
which assignes a particular and unique proper time t, to any given r. If the medium is time- 
invariant, then t(r) = t, /i(r) = /io, e(r) = Eq, reducing (18) to the well-known free-space 
radiation formula (see e.g. [7]): 

x'GR 3 to 

+ / 7 | 1 ^ [3Q-I]^^dx- (20) 
J 47r|x-x'| 2 c £ 

x'GR 3 

y 47r|x-x'| cfeo 

x'GR 3 

and the retarded time is t T = t — |x — x'|/co. Comparing (20) and (18) we see that the 
difference between the time-varying and time-invariant (constant-impedance) background 
media is in the presence of the time-weghting and in the modulation of the source current 
density. It is interesting to note that although the formally defined velocity of light in a time- 



varying medium would be also time-varying, i.e. c = 1/ 'y/e(t)fi(t) = c /a(t), the radiation 
formula (18) features the time delay factor |x — x'|/co, see (19), containing the usual constant 
speed of light in vacuum. In other words, the propagation of light in a time- varying constant- 
impedance medium may be viewed as a propagation in vacuum where the causal time-delay 
is measured in terms of the "universal" r-time. The actual measurements however are done 
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in terms of the usual t-time. Hence, to translate the abstract r-domain picture into the t- 
domain reality, we need to find the inverse mapping t(r) and use it in the given functions a(t) 
and J(x, t) thus obtaining the functions a(r) and J(x, r) appearing in the radiation formula 
(18). Then, after performing all the necessary mathematical operations with respect the 
r- variable in the right- hand- side of (18), we shall try to recognize the t-domain solution. 

Although, the inverse function t(r) exists for any a(t') > 0, to < t' < t, it can be 
found explicitly only if the integral in (5) can be evaluated analytically and the subsequent 
algebraic equation can be analytically solved for t. Consider a linear time-dependence: 



Then, from (5) we find 



a(t) = a + bt. (21) 



Qp + bt bT a , , 

t{T) = e - —, (22) 



and 



a(r) = (ao + bt )e bT . (23) 

A time-harmonic source J(x, t) = I(x) cos(ujt) will have the following representation in the 
r-domain 



J(x, r) = I(x) cos(a;t(r)) = I(x) cos 



+ bto b ao 



(24) 



In the far-field zone the electric field strength is proportional to the last term in (18). 
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Explicitly we obtain 



E(x,r) 



'eR 3 



x'GM 3 



+ 



x'eK 3 



/ 



x'GR 3 



+ 



/ 

'eR 3 



I] I(x')9 r {(a + bt )e bT cos [u (^e bT - f)] } t= ^ , 

(J..X. 



bT 



47TCo5 |x - x'|(a + bt )e' 



- I] I(x')fee 6rr cos [uj ("°±Mo e &r r - f)] , 



bT 



47rcQe |x — x'|e' 



- I] I(x>(a + 6t )e 2fcTr sin [w (M+^i e ^ _ «a)] ; 

(J..X. 



47TCo£o|x - x'|e bT 
— I] I(x')6e ^"' x ~ x 'cos 



47rco£o|x — x'| 



- I] I(x'V(a + bt )e bT e 2 4 |x ~ x ' 1 sin 



-^iidx' 



.,; I «2±Wo e 6T e ^ |x x '' 



47TCo£o|x — X'| 



^ildx' 



(25) 



Where, using (22) and (23), we recognize the t-domain result as 



E(x,t) 



- I] I(x>e~^ |x ~ x '' cos (w't + ^(u/ - u)) , 

— dx 



47TCq£o|x — X'| 



+ 



/ 



- 1] l(x!)uj'a(t)e 



-|x-x'| 



4:7VCqE \x — x' 



sin + ^ 



(26) 



In this formula we have a new frequency 



/ --Mx-x'| 

uj = oje c o 



(27) 



showing that one of the effects due to the time- varying background is a shift of the observed 
frequency depending on the distance between the source and the receiver. In particular, a 
medium varying linearly with time at a rate b introduces a frequency shift varying exponen- 
tially with distance (27) - blueshift for negative b, red-shift for positive. We also notice an 
additional exponential distance factor in (26) - decay for redshift, amplification for blueshift. 
Tracing back to (18) we see that the latter parametric amplification is due to the modulation 
of the source by the a(r)-function, which produces other interesting phenomena as well. For 
example, if the time- variation of the current is such that 



d T [a(r)J(x, t)] = 0, i.e. a(r)J(x, r) = const, 



(28) 
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then the far-field term in (18) disappears, meaning that the field will die-off very rapidly 
with the distance from the source, and that the time- varying current, which would normally 
radiate, may be effectively silenced by a suitably time-varying homogeneous and lossless 
medium. 

The other related but opposite phenomenon is the radiation by a DC-like current. Nor- 
mally, a constant current does not radiate (into the far- field zone), neither does a charged 
particle moving at a constant velocity. Whereas, in a time-varying medium, a constant cur- 
rent would produce a far-field radiation proportional to the magnitude of the current and 
with the time-characteristics depending on a(t). Indeed, 

d T [a(r)J(x)] = J(x)<9 r a(r) ^ 0. (29) 

Obviously, it is the medium, which is the true cause of radiation in this case. 

All of this confirms that a time-varying medium, even a spatially homogeneous one, 
may produce a significant space-time modulation of the electromagnetic field. Although, a 
medium with constant impedance considered so far is a very rare thing, at least one physically 
interesting instant of such a medium is found in General Relativity Theory (GRT) [8, 9], 
where a linearly expanding universe corresponds to the linear a(t) discussed above (apart 
for the current-modulating factor, which is a 4 (r) in GRT). In particular, application of the 
present analysis to the problem of Hubble's redshift law shows that a linearly expanding 
universe produces an exponential redshift, not a linear one as it is currently believed. In 
general, it may be shown, using only (22) and (19), that any two temporal features of the 
time-derivative of the modulated current density, say the zeros or extrema of the signal, 
radiated by the source at times tf° and t| rc , with an interval At src = tf c — tf c between them, 
will be detected by the receiver at times t\ ec and F 2 ec with a different interval At rcc = t™ c — t| ec , 
namely: 

A -i-rec 

At src ' 

where x src and x rcc are the spatial locations of the source and the receiver, respectively. This 
is simply a generalization of the frequency shift law (27). In combination with the radiation 
formula (18) and the time-delay expression (19) this law provides an intuitively appealing 
picture of the expanding universe. First of all, there is no relative motion of the source and 
receiver, as the position vectors in the final radiation formula (18) are time-independent. 

9 



The speed of light also does not change as far as the r-domain picture is concerned. Yet, 
in the real t-domain world we observe the stretching of the time-interval between certain 
characteristic events in the received signal according to (30). 

Imagine two clocks, which are intially at the same location in space and are synchronized. 
Suppose that each second these clocks radiate a short impulse of light. Now, we place one of 
these clocks some distance away and observe that the interval between the impulses arriving 
from that clock is longer than one second (impulses have become broader in time too). We 
are naturally inclined to think that the remote clock started to walk slower for some reason, 
and decide to go and check it. To our surprise, however, upon our arrival we find that 
everything is fine with this remote clock, and it is the other one which walks slower, judging 
from the impulses that are coming from it. We begin to suspect that there is a redshift 
involved due to the relative motion of the two clocks. Yet, when we put a ruler between 
them and observe the measured distance for some time, nothing seems to happen, and the 
clocks appear to be at rest with respect to eachother. This, in short, is what happens in an 
expanding universe. 

Due to its importance in cosmology, especially for understanding the process of inflation 
proposed within the Big Bang theory, we shall consider here the case of exponentially growing 
a(t) as well. Suppose, 

a(t) = a e bt . (31) 

Then, from (5), the inverse function t(r) is found to be 

t( T ) = -jln \e- bt0 - a br] . (32) 
b 

It is obvious that this inverse mapping is only defined for r < e~ bt °/(a b), since otherwise 
the physical time t would become complex. Therefore, according to (19), we may conclude 
that the light impulse radiated at r will be confined to the spatial radius 

R = *L e -bt _ CqT) (33) 
a b 

effectively excluding any communication between points further than R apart. Alternatively, 
one could speculate that the laws of light propagation are different beyond that radius, which 
is highly unlikely. Computing the ratio of time intervals between the corresponding features 
(e.g. zeros) of the electromagnetic signal at the source and receiver locations we get the 
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following expression 



At rcc 



(34) 



At src 



where 



r 



s 



e- bt " - a 6r| rc 

e- bt " - a brf TC ' 

e- bt " -a 6(r| rc + |x-x'|/c ) 

e- w " - a & (7f c + |x - x'I/cq) ' 



(35) 



This proves wrong the common assumption that an exponentially expanding universe in- 
troduces an exponential cosmological redshift. The redshift appears also to depend upon 
the interval between Tf rc and t^, i.e. is inhomogeneous - different frequencies will have 
different shifts. 

III. TIME- VARYING IMPEDANCE 

The problem for a general time-varying impedance with asynchronously varying permit- 
tivity and permeability so far resists analytical treatment. However, for a class of time- 
varying impedances, where /i(t) is a constant, and e(t) is a fourth-order polynomial in time, 
a useful analytical result can still be obtained. The solution can even be put in a form of a 
radiation formula similar to (20) and (18) revealing important differences with the stationary 
medium and the constant-impedance cases. The technique is due to A. B. Shvartzburg, who 
considered mainly the source-free one-dimensional Maxwell's equations. Here this method 
is extended to arbitrary pulsed sources in three dimensions. We denote: 



e(t) = e a(t) = e u 2 (t), ji{t) = /i = const. 



(36) 



We again introduce a new (slightly different) time-like variable: 




(37) 



so that the partial derivatives are related as 



d 1 d 
dt ~ ^dr' 



(38) 
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In addition to this variable change, a new unknown function is introduced as 

F(x,t) = ^. (39) 

To deal with the three-dimensional vector case some preliminary work needs to be done, 
consiting of deriving the (x, r)-domain compatibility relation for the F-field. It can be 
obtained from the (x, r)-domain Maxwell's equations, and turns out to be 



V-F(x,r) = -_ -L=(V.) / u(r')J(x,T')dr'. (40) 

y/u[T) J 


Now we go back to the (x, t)-domain Maxwell's equations and eliminate H, thus obtaining 
the second order vector wave equation for D: 

v 2 (t) v 2 (t) 
V x V x D(x, t) + ^<9 2 D(x, t) = ^<9 t J(x, t). (41) 

Performing both the function and the variable change, and using (40) we arrive at 

F(x,r) 



AF(x,r)-i«9 T 2 F(x,r)-i 



\u{t)dlu{t) - \ {d t u{t)f 



^4^9 r J(x,r)-^==V(V.) / u(r')3(^r')dr'. 



(42) 



Now, and this is the essence of the Shwartzburg method, we shall limit the time-dependence 
of the medium, i.e. the class of functions u(t), to those satisfying the nonlinear equation 

\u{t)dlu{t) - \ {d t u(t)f = -L (43) 

where T is a real constant with the physical dimension of time. In that case (41) reduces to 

AF(x,r)-^F(x,r)-^F(x 1 r) 



m 1 r ^ 

= ^fid T J(x, r) - ^=V(V- ) / u{t>) J(x, t') dr', 

v Q 

which can be solved analytically for any source. The functions u(t) leading to this simplified 
equation are the solutions of (43) and have the form: 

u(t) = a + bt + ct 2 . (45) 
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This polynomial function allows to model both rising and decaying profiles of e(t) = e u 2 (t), 
while parameter T determines the temporal scale of variations. 

The solution of (44) is obtained via the same path as in the previous section. First, we 
transform the equation to the (k, s)-domain and solve it there. Then, we carry out the 
inverse Fourier transform and obtain the following (x, s)-domain solution: 

e -7|x-x'| 



F(x, S )= J J^__*, |3 [SQ-lWx'^dx' 

x'eR 3 

+ / 4llx-xf [3Q-lJi(*'.«)d* / 

x'eR 3 

/ 2 - 7 |x-x'| ^ 
f- -QJ 1 (x', s )dx / 
47r|x — x'| 

x'eR 3 

r e -7|x-x'| i a 

~ / ~a — i T\— J2(x',s)dx / , 

J 47r|x-x / |co 

x'eR 3 

where the currents are the (x, s)-domain images of 

T 

J!(x,r)=^i— / M (r')J(x,r')dT', 



(46) 



o 

J 2 (x,r) =a/m(t)<9 t J(x,t). 



(47) 



The difference with (17) is mainly in the mathematical form of the parameter 7, which is 
now 

7 = -v^ + l/T 2 . (48) 
Co 

Obviously, for 1/T 2 = 0, corresponding to the infinitely slow variation of e(t), we recover 
the free-space case 7 = s/c . To transform the result to the (x, r)-domain, we recall the 
following known Laplace transform pairs: 

e~ k ^+^, (k>0) -> 6(t — k) — — — — J 1 (aVW^¥)H(t - k), (49) 

V£ 2 — & 2 



, (k>0) -> ^(aVt 2 - k 2 )H(t-k), (50) 



Vs 2 + a 2 

where i/(t — fc) is a unit step function located at £ = k. Using these transforms we can write 
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the result as 



F(X ' T)= / 4n\x-xV W-'l-M*'-^*' 

^ / i ^bf[3Q-l]/ 9l (x-x', T -r')J 1 (x',/)dr'dx' 



where 



x'eM 3 



'gr 3 



x'ei 



+ 



Cq y 47r|x — x'| 

x'ei 3 



d 2 T + 





1 

Ji2 



Ji(x',t') dr' dx' 



Ji(x', r r ) dx' 



i_ J 1-Q j gi (x-x', T -T') 

-J 2 (x / ,r r )dx / 



'eR 3 
1 f 1 



Ji(x , ,r')dr / dx' 



Cq y 47r|x — X 



+ cfr / ^/^ 1 ( x ~ x '' r ~ r/ ) j2 ( x/ ' r/ ) dr/dx/ ' 



'eR 3 o 



(51) 



5i(x,r) 



-2 _ |x| 2 



C 



(52) 



and the retarded time is again given by (19). Notice that this retarded time appears both 
explicitly in some terms in (51) as well as implicitly in the Green's functions as a unit step 
multiplier, thus the causality in the (x, r)-domain is fully preserved. Computing out the 
time derivatives 



d 2 + — 



T 

Ji(x,r) =v(t) J M(r , )J(x,r , )dr' + w 1/2 9 r J(x,r), 



v(t) 



o 

1 3d T u dlu 



(53) 



T 2 u l/2 4m 5/2 2M 3 / 2 ' 
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we can explicitly show all the terms in (51) as 
E(x,r) = 

/ 



>m 3 o o 

t' 



i 



3Q 


— I 


1 1 


An x 


-x'| 3 


e(r)J 
o 


3Q 


— I 


1 


47T X 


-x'| 2 


Cqe{t) 


3Q 


-I 


1 


47T X 


-x'| 2 


c e(r) 


3Q 


— I 


1 


An x 

c 


-x'| 2 

1 


c £(t) 
1 



ffl (x- X ',r-r') / ^ I7 ^A 7 ^J( x ',r'')dr''dr'dx' 



+ / , , Q ^ 7 , / £o(x - x', r - rO / ^V2 (r) ^ (r>(r / J(x / ^ dr » d/ ^ 
J 47r|x - x'| 2 c £(r) J J 



'GK» 

T 



./ Ltt|x - x'|- <- () :( t ) 





(54) 



47r|x — x'| CqE{t) 



I -^-r^ j Ml/2 ( r ) w ( r r)«(r / )J(x / ,r / )dr / dx / 

o 

t' 

,V2 
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There are some similarities here with the constant-impedance case considered in the previous 
section. First of all we see both the time- weighting and the modulation of the source. 
However, there are also obvious differences. The time weighting (37) contains now the square 
root of a(t). Hence to achive the exponential spectral shift (27) we need u(t) = a + bt, i.e. 
a quadratic variation of e(t): 

e(t) = e (a + bt) 2 . (55) 

Further, there are new terms containing r-convolutions with the Green's functions (52). 
These terms represent causal (since r-integration starts only after the wavefront has reached 
the observation point) echoes due to the scattering of the field on the temporal variation 
of the medium impedance. Finally, we notice some additional longitudinal contributions in 
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the far-field zone - last four terms in (54). Normally, the far-field contribution is purely 
transverse, i.e. contains only one term with Q — I, see (20). Whereas, here we have three 
additional terms with Q decaying at the typical far-field |x — x'| _1 -rate. Of course, two of 
them are medium- induced echoes. However, one is definitely produced by the source alone. 
This means that the time-varying background could be used to propagate the longitudinal 
near-field features of the source into the far-field. 

IV. CONCLUSIONS 

We considered the effects of two types of homogeneous time-varying background media 
on the electromagnetic field radiated by an arbitrary source. The first type is a medium with 
synchronously varying permittivity and permeability, where the impedance stays constant at 
all times. Such a medium describes the time-vaying spatially flat gravitational background, 
e.g. the expanding universe. The obtained analytical solution shows that under a change of 
time-variable, the problem is similar to the radiation in a stationary background medium. 
Yet, the new time- variable and the modulation of the source term produce significant changes 
in the space-time evolution of the field. We have shown that a linear variation of the medium 
with time produces a "cosmo logical" frequency shift, which is an exponential function of the 
distance between the source and the observer as well as an additional exponential factor 
affecting the amplitude of the received signal. We have also derived the "logarithmic" 
inhomogeneous frequency shift for a universe expanding exponentially in time. 

The more down-to-earth case of the constant permeability and time-varying pemittivity 
appeared, in fact, to be more difficult and was treated using the Shwartzburg method, which 
involves the change of the time- variable as well as the change of the unknown. This method 
gives an analytical solution for a limited but useful class of time-varying permittivities. The 
obtained solution shows, in particular, that the medium with a quadratic temporal change 
in permittivity would also produce an exponential distance-dependent shift in the observed 
frequency. In addition, the received signal is shown to contain causal echoes induced by the 
scattering of waves on the time-variations of the impedance. The modulation of the source, 
which occurs in both constant and time-varying impedance cases may be used either for 
silencing or for amplification of certain sources. Surprisingly, we find that in a time-varying 
background, a DC-like constant current radiates, whereas a specific time-varying current 
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does not. 
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